16.451 Lecture 18: More on Beta Decay 


8 Nov. 2005 


Recall from last class: we calculated the partial decay rate for a final state e- 
with momentum p using Fermi's Golden Rule: 
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where G is a generic “weak coupling constant" for the process. 

• If the electron and neutrino couple to S = 0, then it is called a Fermi decay and G - G v . 

• If they couple to S = 1, it is a Gamow-Teller decay and G = G A . 

• For the neutron decay, both 5 = 0 and S = 1 conf igurations are possible, and S = 1 is 3 
times more likely, so in that case G 2 = G v 2 + 3 G A 2 . 


Electron momentum spectrum: 


N(p)dp 


N 0 *if 


(const.) x p 2 q 2 dp 


N(p) = (const.) x p 2 q 2 


approx: K R = 0 

/ 

(const) x p 2 (Q-K e ) 2 



1 






Comparison to Reality: (Krone, fig. 9.3: e+ and e- decays of 64 Cu) 3 


Too many low energy e 



Too few low energy e + 


Discrepancy: neglect of Coulomb effects in the final state. 


Key point: Coulomb distortions of the energy spectra occur AFTER the electron/positron 
are emitted in the weak decay process. 

Modified density of electron/positron states: 



"Fermi function", depends on 
the charge Z' of the "daughter 
nucleus" (final state) and the 
electron/positron momentum 


Approximate correction factor for |3± decay: 
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► Modified electron/positron spectrum prediction: 
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Fermi-Kurie Plot 


Idea: for "allowed decays", corresponding to our approximation: e ip R .r/n = l 

inside the nucleus, the electron energy spectrum can be "linearized" if one 
accounts for the Coulomb distortion via the Fermi f unction F(Z',p): 



66 Ga -> 66 Zn decay, Phys. Rev. 129, 1782 
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Neutrino Mass effect: 3 H — > 3 He + e + v e (Q = 18.6 keV) 


Idea: shape of the electron energy spectrum near the endpoint (Q) is sensitive to the 
mass of the electron antineutrino: 
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recall: Q = K R + K e + K y 


When K e = Q, K R = K v ^ 0. if m v ^ 0, then 
in this limit, mass effects are most pronounced. 



energy E [keV] E— E a [eV] 


Figure 2: The electron energy spectrum of tritium j3 decay: (a) complete and (b) narrow region 
around endpoint Eq. The (3 spectrum is shown for neutrino masses of 0 and 1 eV. 
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figure from "KATRIN" proposal, 2001: (Karlsruhe Tritium Neutrino expt.) 
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Figure 4: Results of tritium rt decay experiments on the observable over the last decade. 


• Best direct upper limit: m v < 2.2 eV 

• from Sudbury neutrino observatory and other experiments, we have convincing 

indirect evidence of nonzero neutrino mass that is much smaller than this 


Web sites: http://cupp.oulu.f i/neutrino/nd-mass.html , http://www-ik.fzk.de/~katrin/ 


Systematic test: total rate for beta decay 


Our formalism determines X ift which is the rate (s 1 ) to a particular final state electron 
(or positron) momentum p: 
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(± refers to ff decay modes) 


The total decay rate is obtained by integrating X if over all allowed e 1 momenta p: 
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"Fermi integral", f (Z',Q) 


Key point: apart from the nuclear matrix element , the variation in decay rates for 
different unstable nuclei should only depend on the Fermi integral , which we can 
calculate independently. 

We can use this to test our weak interaction theory! 


4 





Krone, Figure 9.& Dimensionless Fermi integral 


By convention: f(Z',E 0 ) = -f - ] p\E 0 -K e ) 2 F(Z',p) dp , E 0 

XYYC J 



Comparative half lives: 


10 


By convention, the half -life, t 1/2 =T In 2 , WZth T =1/ A is used as a comparison 


standard for different nuclear beta decays: 

we had: X (s 

_1 > = , £, 3 >< Kdf f(Z',Q) 

2n n c 

rearranging, we get: 

dimensionless integralx^^ 
from last slide 
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Notice: The only difference in the "ft" value between different nuclear beta 
decays is the value of the nuclear matrix element. 



If | M nuc | | 2 = 1 ("superallowed" case in nuclei), the ft values can 
be used to determine the weak coupling constants G - (G v , G A ) 


Special case: "superallowed" decays in nuclei with initial and final nuclear states 
0+ -» 0+, e.g. 
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O -> X yN + e + + v e 


must have S = 0 for the leptons -> pure Fermi decay... 
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Superallowed Fermi 0+ -» 0+ decays: world's best data for light nuclei 
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• all have the same ft value ~ 3100 sec 

• determines the weak coupling constant for Fermi decays: 

G v = (1.1358 ± 0.0004) x 10~ 5 ( he) 3 / GeV 2 
(And G a /G v = -l 25, more la ter. . . .) 


Can we understand beta decay rates in general? 


first page of Krane, Appendix C: (symbol s stands for electron capture/ p+ decay) 


27 isotopes: 8 p - decays, 6 p + decays, spanning 16 orders of magnitude in rate! 
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Some anomalies: 


1. According to our theory, the very slow decay: (1.6 x 10 6 yrs) 


^Be (0 + ) — » ^B(3 + )+ e“ + 


should not occur at all, because angular momentum does not add up, i.e.: 

6 * 3 + (6 or 1) 

2. Another example: (16.1 hr) 


35^r(l ) — > 34 Se ( 0 + ) + e + + v e 


This should not occur because the wavefunctions in the nuclear matrix element 
have opposite parity, so the integrand is odd and should vanish: 


M nuclear = JVf(0 ^,'(0 d 3 r = 0 ??? 


Forbidden Decays: 


These are two examples of forbidden decays - they cannot proceed under the allowed 
approximation, since 

M if = G!v'* P , f (r)fe(r)f v (r)V'n,i(> : ) d 3 r = 0 if f e (r) £{?) = L 
Is there some other way they can occur? 

Reconsider the electron and antineutrino wave function as a multipole expansion: 

v f e {r) f v (r) = e l P R ' r ' n = fi L (2L+l) j L (p R r/h ) P L ( cos0) 

L = 0 


P^r 

Pr ^•^ipr+rnn Jl = spherical Bessel Function 

' Ox 

recoil P L ( c ° s 6) = Legendre polynomial 

Q neutrino 
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Beyond the "allowed approximation": 
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ip R >r/n _ ^,-1(21+1) i L (p R rlh) P L (cos#) 


spherical Bessel functions: 

• / \ sinx . . . sinx cosx 

Jo( x )= — ; Ji( x ) = — 2 •••• 

x X X 

with x = p R r/h 

for successively larger L, they become 
more significant for larger recoil momentum 
-> this will change the momentum 
dependence of our prediction ! 

Legendre polynomials: 

P 0 =l, P 1 = x, P 2 =-(3x 2 -1).... 
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with x = cos& 

these introduce a new angular dependence 


to the integrand for M if equivalent to angular momentum L 


Effect on beta transitions: 
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• angular momentum coupling for the multipole order L, together with S and nuclear 

angular momentum allows previously impossible reactions to proceed 

• multipole term has parity (-1) L , which allows nuclear states of opposite parity to 

be "connected" by the beta decay operator 

• momentum dependence of the matrix element varies as (Pj?r / ^) L ... 

since this is small, the lowest multipole order L that satisfies the conservation 
laws will dominate the transition 


rate- | M | 2 - (p R r/h) 2L = (0.01) 2L -> dramatically smaller for large L 

momentum dependence also affects the shape of the spectrum; Kurie plots 
are not linear unless "shape factors" are taken into account.... 


• naming convention: 
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L = 3 
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Class if i cat ion : all known decays.... 


Allowed transitions 


£ 
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variation within a 
class is due to the 
nuclear matrix element 


Observed log 10 ft value 


Figure 5-8: Systematic* of observed log ft values. The grey area in the upper 
panel shows 718 cases of 0 + *=* 1 + allowed transitions, and the remaining 1741 
cases of other allowed decays are shown by the white histogram. The peak of 
the distribution for the 24 cases of 0 + -» 0+ superallowed decay is indicated 
by the arrow. 


Classif ication of decays: general 


Nuclearcase: ^ Z' Y + e + V e 
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(S, L) 


Conservation laws: j _ J ^ + S + L 

n i = 7i f (-1) L 

with S = 0 (Fermi) or S = 1 (Gamow-Teller) 

Smallest value of L that is consistent with conservation laws will 
dominate the transition. 



Example: p+ decay of 18 Ne 



Branching ratio (BR): the f raction of decays that go to a particular final state. 
In this case, 7 tota) = 1/x = 0.667 sec 1 ; 7 = + 7 2 + 7 3 , with 7; - BR(i) X tota) 

Transition 1: 0 + ->■ 0" This is a first forbidden GT decay, with the slowest partial rate: 

6 = 6 + S+L) (+) = (-)x(-l) L L = 1, S = 1 
Transition 2: 0 + 0 + This is an allowed Fermi decay. 

0 = 0 + S+L; (+) = (+)x(-1) l -> L = 0, S = 0 

Transition 3: 0 + -»1 + This is an allowed Gamow-Teller decay 

6 = 1 + S + L ; (+) = (+)x(-l) L 1 = 0, S=1 




